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ABSTRACT: The present paper discusses solution of dispersion phenomenon by using Homotopy analysis method. 
Solution represents concentrations of any contaminated or salt water disperse in homogenous porous media saturated with 
fresh water. The solution of non-linear partial differential equation has been obtained in term of series solution of 
exponential function ofX and time T under assumption of guess value of concentration of contaminated or salt water. Here 
solution converges for parameter s = 1 for embedded parameter ti = 0. 1 . The graphical presentation is given by using 
Maple coding. It is concluded that concentration of contaminated or salt water dispersion is decreasing when distance X as 
well as time T increasing and convergence of the solution has been discussed. 
Keywords: Advection, Diffusion, Homotopy analysis Method 



I. INTRODUCTION 

The present paper discusses the solution of longitudinal dispersion phenomenon, which arising in the miscible fluid 
flow through homogenous porous media. The problem of solute dispersion during underground water movement has 
attracted interest from the early days of this century [1], but it was only since 1905 in general topic of hydrodynamic 
dispersion or miscible displacement becomes one of the more systematic studies. The phenomenon of the dispersion has 
been receiving good attention from hydrologist, agriculture, environmental, mathematicians, chemical engineering and soil 
scientists. The specific problem of fluid mixing in fixed bed reactors has been investigated by Bernard and Wilhelm [2]. 
Kovo [3] has worked with the parameter to be modeled in the longitudinal or axial dispersion coefficient D in chemical 
reactors model. 

The fundamental interest of this paper is to find concentration of contaminated or salt water. The term concentration 
expresses a measure of the amount of a substance within a mixture. The dispersion process is associated with molecular 
diffusion and mechanical dispersion. Molecular diffusion is the spreading caused by the random molecular motion and 
collisions of the particles themselves and mechanical dispersion is the spreading of a dissolved component in the water phase 
by variations in the water velocity (i.e. flow of a fluid). These two basic mechanisms molecular diffusion and mechanical 
dispersion cause a concentration front of fluid particles to spread as it advances through the porous media. These two 
combine processes of molecular diffusion and mechanical dispersion are known as hydrodynamic dispersion or dispersion. 
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Fig 1: The geometry of microscopic pores, where velocity distributions in different pore size. 



When groundwater flows, the actual microscopic velocity in the pores varies widely in space even when the Darcy 
macroscopic velocity is constant. The result is more intense mixing, which is called hydrodynamic dispersion. Figure 1 gives 
a schematic view of the trace movement on macroscopic level. This phenomenon can be observed in coastal areas, where the 
fresh waterbeds are gradually displaced by seawater. This phenomenon plays an important role in the seawater intrusion into 
reservoir at river mouths and in the underground recharge of groundwater. 

Most of the works reveal common assumption of homogenous porous media with constant porosity, steady seepage 
flow velocity and constant dispersion coefficient. For such assumption Ebach and White [4] studied the longitudinal 
dispersion problem for an input concentration that varies periodically with time. Hunt [5] applied the perturbation method to 
longitudinal and lateral dispersion in no uniform seepage flow through heterogeneous aquifers. Mehta and Patel [6] applied 
Hope-Cole transformation to unsteady flow against dispersion of miscible fluid flow through porous media. Marino [7] 
considered the input concentration varying exponentially with time. Eneman et al.[8] provided analysis for the systems 
where fresh water is overlain by water with a higher density in coastal delta areas. Experimental evidence of such lenses 
was given by Lebbe et al. [9], and Vandenbohede et al. [10]. Meher and Mehta [11, 12] studied the Dispersion of Miscible 
fluid in semi infinite porous media with unsteady velocity distribution using Adomain decomposing method. 

The present paper discusses the approximate analytical solution of the nonlinear differential equation for 
longitudinal dispersion phenomenon which takes places when miscible fluids (contaminated or salt water) mix in the 
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direction of flow. The mathematical formulation of the problem yields a non linear partial differential equation. The 
analytical solution has been obtained by using homotopy analysis method. 

II. Statement of the problem 

Considering dispersion of contaminated or salt water with concentration c(x,t) flowing in x-direction, dispersion 

taking place in porous media saturated with fresh water. Hence it will be miscible fluid flow through homogenous porous 
media. Therefore, it will obey the Darcy's law, which dates back to 1856 [13]. The following assumptions have been made 
for present analysis (Schidegger 1954, Day 1956, deJony 1958) [14,15, 16]: 

• The medium is homogenous. 

• There is no mass transfer happen between the solid and liquid phases. 

• The solute transport across any fixed plane, due to microscopic velocity variation in the flow tube, may be 
quantitatively expressed as the product of a dispersion coefficient and the concentration gradient. 

To find concentration of the dispersing contaminated or salt water as a function of time t and distance x, as the two 
miscible fluids flow through homogenous porous media. Since the mixing (contaminated or salt water and fresh water) takes 
place both longitudinally and transversely. Dispersion adds a spreading effect to the diffusion effects. Science dispersion is 
driven by the mean flow of the water, the dispersion coefficients related to the characteristic length or pore length L. In three 
dimensions, the spreading caused by dispersion is greater in the direction of the flow than in the transverse direction. One 
dimensional treatment of these systems avoids treatment of a radial or transverse component of dispersion. We only consider 
the dispersion phenomenon in the direction of flow (i.e. longitudinal dispersion), which takes places when miscible fluids 
flow in homogenous porous media. 

III. Mathematical Structure 

According to Darcy's law, the equation of continuity for the mixture, in the case of incompressible fluids is given by Bear 
[1]. 

f + v.H=o (i) 

Where, p is the density for mixture [ML" 3 ], t is time [T] and V is the pore seepage velocity vector [LT" 1 ]. 

The equation of diffusion for a fluid flow through a homogeneous porous medium, without increasing or 
decreasing the dispersion of contaminated or salt water is given by 



dt v i 



pD 



3, 



(2) 



Where C is the concentration of dispersing contaminated or salt water per unit volume [ML" ] and D is the tensor 
coefficient of dispersion with nine components [L 2 T -1 ]. 

Since contaminated or salt water is flowing through a homogeneous porous medium at constant temperature, p may 
be considered as constant. Then 

v-v = o (3) 

where V is velocity of contaminated /salt water dispersion. 
Hence Eq. (2) may be written as 

^ = -v.vc+v-\dvc] ( 4 ) 
dt L J 

When the seepage velocity V of contaminated or salt water is dispersing along the x-axis, then the non-zero components will 
be n _ n _ L , (coefficient of longitudinal dispersion and L is length of dispersion in flow direction) and 

*-() 

D n = D T (coefficient of transverse dispersion) and other /> are zero. From this assumption the equation (4) becomes, 

dC dC „ d 2 C e*\ 
— = -u — + D L — - yj > 

dt dx dx~ 

where U is the component of seepage velocity of contaminated or salt water in x direction which is function of x and t 
and d l > 0 ■ It has been observed that seepage velocity U is related with concentration of contaminated or salt water 
dispersion. We assume that seepage velocity u is directly proportional to c(x,t) U 1] 

tt - C M (6) 
Q 

where \lr is constant of proportionality and the guess approximation of the concentration of contaminated or salt water 
dispersion. To get dimensionless form of equation (5) using the dimensionless variables 

X=^,andT = -t 
L L 
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then equation (5) can be written as, 

^ = -C — + f — ' where e = DlC ° and e efO,ll,0<X <1,0<7/<1 (7) 

8T 8X dX 2 L 

Since concentration C is decreasing as distance X increase for T > 0. It appropriate to choose guess value of concentration 
solution as, [12] 

U (X,T;e) = ^\-^XTje- x +e m ( 8 ) 

Hence, the equation (7) together with boundary condition (8) represents the governing non-linear partial differential 
equation for concentration of the longitudinal dispersing material of miscible fluids flowing through a homogenous porous 
medium. 

IV. The solution with Homotopy Analysis Method 

For one dimensional non-linear partial differential equation for longitudinal dispersion phenomenon, we assumed that the 
concentration c[X,T) °f tne dispersing contaminated or salt water, at time T=0 is expressed as, 

□ ( X T s = 0) = f 1 - — TX J e~ x + s'" > wnere = 0 for the concentration of contaminated or salt water (9) 

Now we apply the homotopy analysis method into the longitudinal dispersion phenomenon during miscible fluid flow 
through homogenous porous media. We consider the equation (7) as nonlinear partial differential equation as 

□ [□(X,T;ff)] = 0 (10) 
Where □ is a non-linear operator, □ (x,T;e)is considered as unknown function which represent the concentration C of the 
dispersing contaminated or salt water at any distance X for given time T > 0, for 0 < £ < 1 . We use auxiliary linear 
operator jjj T-ejJ = — (^'^' g ) and initial approximation of concentration of dispersing contaminated or salt 

water ^ ^ r) = |l— — TXJe ~ x t0 constmct the corresponding zero* order deformation equation. As the auxiliary linear 

operator 3 which satisfies 3[c 4 ] = 0> where c 4 is arbitrary constant. This provides a fundamental rule to direct the choice of 
the auxiliary function H(X, T) the initial approximation c o (X,T), an d tne auxiliary linear operator 3, called the rule of 
solution expression. Establish the zero-order deformation equation of longitudinal dispersion phenomenon as [20], 
(l-s)3[D (X,T;s)-C 0 (X,T)\ = sTiH(X,T)U [□ {X,T;ej\ ( n ) 
where c o (X,T) denote an initial guess value of concentration of dispersing contaminated or salt water of the exact solution 
C(X,T) which is our purpose to find it. Since ^0 an auxiliary parameter and h(X,T\ * 0 an auxiliary function such that 
e e[0,l] an embedding parameter. The auxiliary parameter ti is providing a simple way to ensure the convergence of series. 
Thus it renamed ft as convergence control parameter [20]. Let 3 an auxiliary linear operator with the property that, 
3[0 (X,T;e)] = 0 when c(X,T;s) = 0 

When s = 0 , the zero-order deformation equation (11) becomes 

3[a(X,T;e)-C 0 (X,T)] = 0 (12) 
Which gives the first rule of solution expression and according to the initial g uess c (x r) = 7X je~* ' ^ ^ s 

straightforward to choose 

a (X,T;0) = C 0 (X,T) (13) 
When £ = 1 , since h&Q, h{X,T\ ^ 0th e zero-order deformation equation (7) is equivalent to 

□ [□(X,2V)] = 0 (14) 
which is exactly the same as the original equation (10) provided 

□ (X,T;l) = C(X,T) (15) 
According to (13) and (15) as the embedding parameter £ increases from 0 to 1, solution □ <X,T;s) varies continuously 
from the initial guess value of the concentration c o (X,T) of dispersing contaminated or salt water to the solution 
C(X,T) an d its solution is assumed by expanding □ (x,T;s) m Taylor series with respect to s as, 

U (X,T;s)=D (X,T;0) + '£C m (X,T)e"' ( 16 ) 



Where, Q (x T) _ 1 d'V{X,T;e) 
'"^ ' ' m! ds m 



(17) 
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i.e. the concentration of dispersing contaminated or salt water is function of distance X and time T for any parametric value 
s is expressed as, the concentration of dispersing contaminated or salt water at time T = 0, C 0 (X,T) an d sum of 

concentration of dispersing contaminated or salt water c t (X,T), C 2 (X,T),--- ai different time T for different value of 
parameter s- Here, the series (16) is called homotopy-series; the series (16) is called homotopy series solution of 
□ [□ (X,r;f)] = O anc l C m (X,T) is called the mth-order derivative ofD . Auxiliary parameter ti in homotopy-series (16) can 

be regard as iteration factor and is widely used in numerical computations. It is well known that the properly chosen iteration 
factor can ensure the convergence of homotopy series (16) is depending upon the value of %, one can ensure that convergent 
of homotopy series, solution simply by means of choosing a proper value of h as shown by Liao [20]. If the auxiliary linear 
operator, the initial guess, the auxiliary parameters, the auxiliary function H[X,T) al ' e so properly chosen, the series (16) 

converges at e = 1 ■ 

Hence the concentration of dispersing water can be expressed as, 
C(X,T) = C 0 (X,T) + f j C m (X,T) 

m=l 

And c (X,T) can be calculated by equation (23). This must be one of solution of original non-linear partial differential 

equation (7) of the concentration of dispersing contaminated or salt water problem in homogenous porous medium. 
According to the definition (17), the governing equation can be deduced from the zero-order deformation equation (11), 
define the vector 

C m ={C 0 (X,T),C 1 (X,T),...C n (X,T)} 

Differentiating equation (11) m-times with respect to the embedding parameter £ and then setting £" = 0and finally 
dividing them by m ! , we have the so-called m th order deformation equation of the concentration c(X,T) will be as, [20] 

3[C m (X,T)-%„ l C m _ l (X,T)] = £fiH(X,T)R m [C m _ l ,X,T) ( 19 ) 

1 gT-'O [D (X,T;e)] 



Where 



R m {C m _ v X,T\- 



-1)1 



de'" 



And 



Xm 



(20) 
(21) 



(0, m<\ 
[l, m > 1 

It should be emphasized that c, (X,r)f° r W > 1, is governed by the linear equation (20) with the linear boundary condition 

that came from original problem, which can solved by symbolic computation software Maple as bellow. The rule of solution 
expression as given by equation (8) and equation (11), the auxiliary function independent of £ can be chosen as u (X,T) = 1 , 

[20] 

According to (15) and taking inverse of equation (19) the equation (20) becomes, 

c m (x,t) = z „c,„_, (x,T)+m-' [r„, (c,„_„x,r)] 
R (c „x,r)= i (*■*■=«)] 

"'^ ™" 1 ' (m-l)\ ds'" 
In this way, we get c (X,T) for m=l, 2, 3, ... successively by using Maple software as, 



(22) 
(23) 



C l (X,T) = ^hT(-T 2 X 2 +T 2 X +3TXe x -6Te x +6TX -3T -6Xe x -Yl-Yle x }i 



C 2 (X,T) = Th 



-0.5hXe 2x -0.25TXe 2x +0.052/iT 3 +0.U6hT 3 X V +0.58AT 2 -0.92hT 2 Xe x +hTXe x 
hl.58hT 2 e x +2.5hTe x -he 2x -0.5Te 2x -0.5Xe 2x -he" -0.5hTe 2x +0.19hT 1 e x 



-0.0$3T 2 X 2 e x +0.5TXe x +0.083T 2 Xe x 



-0.0S3hT 2 Xe 2 



-0.25hT 2 X 2 e x +L5hT 



-hT 2 X +0.25hT 3 X 2 -0.32hT 3 X +0.05hT 4 X 2 -Omih^X -0.025hT A X 3 
-0A3ShT'Xe x -e x +0.5hTXe 2x -Q.25Te x +0.33hT 2 e 2x 



(24) 



(25) 



Using initial guess value of concentration from equation (8) and successive 



C(X,T): 



1 



1 — XT \e~ x + — hT 
12 



1 f-T 2 X 2 +T 2 X+3TXe x -6Te x 



,+6TX-3T-6Xe x -\2-\2e x 



+Th 



f -Q5hXe 2x -0.25TXe 2x +0.052hT 3 +0.\46hT 3 X 2 e x +0.58/iT 2 
-0.92hT 2 Xe x +hTXe x +1.5ShT 2 e x +2.5hTe x -he 2x -0.5Te 2x 
-0.5Xe 2x -he x -Q.5hTe 2x +0.19ftTV -0.0837/ 2 X V +Q.5TXe x 
+0m3T 2 Xe x -e 2x -0M3hT 2 Xe 2x -0.25hT 2 X 2 e x +1.5hT 
-hT 2 X+0.25hT 3 X 2 -0.32hT i X +0.05hT*X 2 -O.OllhT^X 
-0.025/;r 4 A' 3 -0.438/tr 3 Xe x -e x +0.5hTXe 2x -0.25Te x +0.33WV 



(26) 
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V. Numerical and Graphical solution 

Numerical and graphical presentations of equation (26) have been obtained by using Maple coding. Fig 2 represents 
the graphs of concentration c(X,T) vs. distance X, for T= 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0 is fixed, and Table I 

indicates the numerical values. The fig 2 and the table 1, indicate the graphical representations of the longitudinal dispersion 
phenomenon of the concentration. The convergence of the homotopy series (16) is dependent upon the value of 
convergence-parameter h [17, 18, 19, 20]. Therefore we choose proper value of the convergence -parameter h = O.lto obtain 
convergent homotopy-series solution [20]. 



C( X T) 

Table I: Concentration of the Contaminated or Salt Water ^ ' ' For Different Distance X For Fixed Time T : 

0.3,0.4, 0.5,0.6, 0.7, 0.8,0.9, 1.0. 



0.1,0.2, 



Distan 


C(X,T) 


C(X,T) 


c(x,t) 


C(X,T) 


C(X,T) 


c(x,t) 


C(X,T) 


C(X,T) 


c(x,t) 


C(X,T) 


ceX 


T= 0.1 


T= 0.2 


T= 0.3 


T= 0.4 


T= 0.5 


T= 0.6 


T= 0.7 


T= 0.8 


T= 0.9 


T=1.0 


0.1 


0.8674 


0.8287 


0.7887 


0.7476 


0.7054 


0.6623 


0.6183 


0.5736 


0.5282 


0.4822 


0.2 


0.7814 


0.743 


0.7037 


0.6635 


0.6225 


0.5807 


0.5383 


0.4953 


0.4517 


0.4077 


0.3 


0.7038 


0.666 


0.6274 


0.5881 


0.5482 


0.5077 


0.4667 


0.4253 


0.3834 


0.3412 


0.4 


0.6338 


0.5966 


0.5589 


0.5206 


0.4818 


0.4426 


0.403 


0.3631 


0.3228 


0.2822 


0.5 


0.5707 


0.5343 


0.4974 


0.4602 


0.4226 


0.3847 


0.3465 


0.308 


0.2692 


0.2302 


0.6 


0.5137 


0.4782 


0.4424 


0.4063 


0.3699 


0.3333 


0.2964 


0.2594 


0.2221 


0.1846 


0.7 


0.4624 


0.4279 


0.3932 


0.3583 


0.3231 


0.2878 


0.2523 


0.2167 


0.1809 


0.1449 


0.8 


0.4162 


0.3828 


0.3492 


0.3155 


0.2816 


0.2476 


0.2135 


0.1793 


0.1449 


0.1105 


0.9 


0.3745 


0.3423 


0.3099 


0.2775 


0.2449 


0.2123 


0.1795 


0.1467 


0.1138 


0.0808 


1.0 


0.3369 


0.3059 


0.2748 


0.2437 


0.2124 


0.1812 


0.1498 


0.1184 


0.0869 


0.0554 
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Fig 2: Represents concentration of contaminated or salt water c(X,T) vs - distance X for auxiliary parameter tl — 0. 1 and 
auxiliary function H(Z,T) = l [20] when T =0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0 is fixed. 



VI. Conclusion and Discussion 

The equation (26) represents concentration of the contaminated or salt water for any distance X and time T >0. It is 
converges for embedding parameter S = 1 and for auxiliary parameter tl =0.1 which is expressed in term of negative 
exponential term of X and time T > 0 . Concentration C will be one from guess value of the exact solution for X = 0 , T= 
0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0. Fig 2 represents the solution for concentration C vs. distance X shows that 
concentration of the contaminated or salt water is decreasing as distance X increasing for T > 0 . From fig. 2 it can conclude 
that for T = 0. 1 concentration of contaminated or salt water is decreasing as distance X increasing and when time is 
increasing and due to different deformation added to C, the concentration of contaminated or salt water is successively 
decreasing exponentially. Since the equation (7) is diffusion type Burger's equation for longitudinal dispersion phenomenon. 
Hence solution is graphically as well as physically consistent with phenomenon. This is physically fact with phenomenon of 
the longitudinal dispersion of contaminated or salt water in homogenous porous medium. 
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